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ABSTRACT 

OS 

OS . We use recent information on geodetic precession of the binary pulsar B1913+16 along with 

measurements of its orbital parameters and proper motion to derive new constraints on the 
immediate progenitor of this double neutron star system. As part of our analysis we model the 
motion of the binary in the Galaxy after the second supernova explosion, and we derive constraints 
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on the unknown radial velocity. We also obtain limits on the magnitude and direction of the kick 
velocity imparted to the pulsar companion during the second supernova explosion. We consider 
the complete set of possible cases, depending on the kinematic age of the system and the 180° 
ambiguity in the pulsar spin orientation. Most interestingly, we find that the natal kick must have 
been directed almost perpendicular to the spin axis of the neutron star progenitor, independent 
of the specific pre-supernova configuration. Such a tight constraint on the kick direction has 
important implications for the physical mechanism responsible for the kick. 



1. Introduction 

Stellar systems with one or more compact objects can provide us with valuable information regarding 
the collapse of massive (;> 8 — 10 M Q ) stars and the formation of neutron stars (NSs) and black holes. Basic 



issues such as the relation of progenitor mass to the type and mass of the compact object as well as the 



effects of mass ejection or fallback are still poorly understood. Furthermore, despite the accumulation of 



pieces of evidence for kicks imparted to NSs at birth in recent years (see van den Heuvel & van Paradijs 1997; 
Burrows 1998; and references therein) the physical origin of kicks and the factors determining its magnitude 



and direction remain unknown. 



Among the different types of binaries with compact objects, double NS systems offer important ad- 
vantages in studying NS progenitors and kicks, mainly because after the second supernova (SN) explosion 
these systems have not evolved in any way other than through energy and angular momentum loss due to 
gravitational wave emission. Therefore their present observed properties can be related to the conditions 
just prior to the second SN in a straightforward way and can be used to constrain the parameters of the 
pre-SN system and the kick imparted to the second NS. Such constraints on the progenitors of double NS 
systems have the capability of testing formation and evolutionary sequences for binary systems and thus 
provide input to models for the population of Galactic double NS binaries. This is of particular interest, 
since merging double NS systems are currently considered as the most promising sources of gravitational 
waves (Schutz 1986) and a possible origin of 7-ray bursts ( Paczyhski [ 1986). 
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The double NS system containing the radio pulsar PSR B1913+16, which was discovered in 1974 by 
Hulse and Taylor ( Hulse fc Taylor "| l975), proved to be the most exciting test laboratory for relativistic 
theories of gravity. General relativity passed this test with flying colors and provided us with detailed 
information about the system. Continuous timing observations of PSR B1913+16 over a period of twenty 



years (Taylor et al. 1976; Taylor, Fowler, fc McCulloch 1979; Taylor fc Weisberg 1982; Taylor fc Weisberg 



1989; Taylor 1994) allowed a highly accurate determination of the spin parameters, the pulsar position on 
the sky, the five Keplerian parameters of the orbital motion, and three post-Keplerian parameters related 
to effects predicted by general relativity. The measured values for the three post-Keplerian parameters are 
in excellent agreement with general relativity and allow a precise determination of the (observed) masses of 
the pulsar and its companion, M p and M c , the semi-major axis of the (relative) orbit, A, and the inclination 
of the orbital plane with respect to the plane of the sky, i. The proper motion of the pulsar on the sky has 



also been measured with moderate accuracy (Damour & Taylor 1991). In addition to spin, astrometric and 
orbital parameters, timing observations at different frequencies has given a precise value for the dispersion 
measure, DM, of the pulsar, which can be converted into a distance given a model for the distribution of 
free electrons along the line-of-sight. Table |l| lists all the measured and derived parameters for the PSR 
B1913+16 binary system relevant to the present study. 

Reports about changes in the relative amplitudes of the prominent leading and trailing components of 



the integrated profile of the pulsar by Weisberg et al. (1989) and Cordes et al. (1990) have been explained 
as a result of a change in the direction of the pulsar spin caused by geodetic precession. This requires a 
misalignment between the orbital angular momentum and the spin of the pulsar. Based on recent observations 



with the 100-m Effelsberg radio telescope, Kramer (1998) also detected a decrease in the relative separation 
of these components, strengthening the case for geodetic precession. Assuming a circular hollow-cone like 
beam, Kramer (1998) was able to derive a value for the misalignment angle of, i.e. — 22(— 4/ + 2)° or 
= 158(— 2/ + 4)°. More recently Taylor (1999) used pulse structure and polarization data measured with 
the 305m Arecibo radiotelescope and obtained 9 = 14(±2)° or — 166(±2)° as the most likely inclination 
of the pulsar spin axis with a safe upper (lower) limit of 22° (158°). In what follows we adopt this upper 
(lower) limit, since the limits on the pre-SN and kick parameters derived here become tighter for smaller 
(larger) values of 0. 

Given our general understanding of the evolutionary history of double NS systems, the spins of the binary 
components and the orbital angular momentum axis before the second SN were aligned due to mass transfer 
in the binary. The inferred misalignment then between the pulsar spin and orbital angular momentum has 
been caused by the supernova kick given to the second NS and corresponds to the angle between the pre-SN 
and post-SN orbital planes. With this new piece of information about the misalignment angle and the proper 



motion measurements at hand we extend the work by Fryer & Kalogera (1997) and combine all available data 
for the PSR B1913+16 binary system. We analyze the motion of the system in the Galaxy and the effect 
of asymmetric SN explosions on binary orbits and we are able to constrain the properties of its immediate 
progenitor as well as its Galactic velocity and the actual birth kick (magnitude and direction) imparted to 
the pulsar companion during the second SN explosion. The constraints on the kick direction have important 
implications for the physical mechanism responsible for the SN kick. 

In the next section we discuss the position and motion of PSR B1913+16 in the Galactic gravitational 
potential in order to estimate its kinematic age and place of birth. In Section |^ we give relations between the 
present binary system and its immediate progenitor. These relations are applied in the final section followed 
by a discussion. 
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2. Motion in the Galaxy 



The velocity of the center-of-mass of a binary system after a SN explosion is determined by properties of 
the kick imparted to the newly-born NS and the pre-SN binary parameters]^] Here, we estimate the center- 
of-mass velocity, Vcm, of PSR B1913+16 after the second SN explosion by making use of its measured 
proper motion. The proper motion of the system at present is related to a number of different velocity 
components in addition to the NS kick velocity. The post-SN velocity is the vector sum of the circular speed 
due to Galactic rotation at the SN location, of the peculiar velocity with respect to this local velocity, and 
of the velocity imparted to the system due to the NS kick and the mass loss associated with the explosion. 
Depending on the age of the system this velocity sum is further altered by the gravitational acceleration in 
the Galactic potential, so its velocity relative to the present Galactic environment may not be representative 
of the velocity of the system right after the second SN explosion. 

Given that the details of the complete evolutionary history of PSR B1913+16 are unknown, to make 
progress in estimating the post-SN Vcm, we make a few simplifying assumptions, which, however, are consis- 
tent with our qualitative understanding of the formation of double NS systems 



1982; Bhattacharya fc van den Heuvel 1991; Lipunov, Postnov, fe Prokhorov [1996; 



rinivasan & van den Heuvel 



Bagot 



1997; Portegiei 



Zwart fc Yungclson | 1998). The initial progenitor of PSR B1913+16 is thought to be a binary with two 



massive stars. At the time of the first SN explosion the system was a relatively wide binary with a typical 
orbital velocity of a few 10 km s . The corresponding center-of-mass velocity imparted to the system at 
the first explosion is of a similar magnitude (Kalogcra 1996). Therefore, we assume that the progenitor 
was in the Galactic plane at the time of the second SN explosion and that its peculiar velocity was small 
compared to the local Galactic rotational velocity (~ 200 km s -1 ; see also Brandt fe Podsiadlowski | l995). 
The observed pulsar is the older of the two NSs, created in the first SN explosion, and has been recycled 
during an accretion phase. We further consider the characteristic age of the pulsar (r c = P/2P ~ 109 Myr) 



to be an upper limit on the time, Tsn, elapsed since the second SN explosion. In fact, Arzoumanian, Cordes 



fe Wasserman | (1998) derived an upper limit of 80 Myr by estimating the equilibrium period of the pulsar 
during the accretion phase (spin- up limit t su ). 



We take as a model for the Galactic potential the model of Kuijken fe Gilmorc (1989) (see Appendix A) 



and integrate the position and velocity of PSR B1913+16 back in time from the present values to the time of 
the second SN explosion. We introduce a Cartesian coordinate system (X, Y, Z) with the Galactic center at 
its origin, the Z-axis perpendicular to the Galactic plane, and the Sun at (Aq, 0, Z@), where X@ = — 8.5±0.5 
kpc and Z Q ~ 15 pc ( feast & Whitelock [ 1997; |Cohcn [ 1995). The Y -axis of our coordinate system is at a 
Galactic longitude I — 90° and in direction of the Galactic rotation. 

The determination of the pulsar's present location in the Galaxy requires its position on the sky and its 
distance from the solar system. While the Galactic longitude, I, and the Galactic latitude, b, are known with 
good precision, I — 49.968° and b = 2.122°, the only distance indicator available at present is the dispersion 
measure based on a model for the Galactic distribution of free electrons. Damour fe Taylor (1991) give 



d = 8.3 ± 1.4 kpc as a distance to PSR B1913+16. The model of [Taylor fe Cordes | (1993) gives a distance of 
7.1 kpc which lies within the error bars given above. For a distance of 8.3 kpc the present position of PSR 
B1913+16 in the Galaxy Xp, is (—3.16,6.35,0.33) in units of kpc. Hence, at present the pulsar is about 
300 pc above the Galactic plane. 



7 Referring to pre- and post-SN systems in the following, we always mean the system immediately before and just after the 
second SN explosion. 
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The proper motion of PSR B1913+16 in right ascension and declination is known with about 10% 
accuracy ( Damour fc Taylor 1991): 

[i a = a cos <5 = — 3.27 ± 0.35 mas yr _1 , (1) 

US = 6 = +1.04 ±0.42 mas yr" 1 . (2) 

The corresponding proper motion in Galactic longitude and latitude is |^] 

Hi = lcosb= -0.60 ± 0.45 mas yr" 1 , (3) 
^ b = b = +3.38 ±0.31 mas yr" 1 . (4) 

The radial velocity of the binary, V r , with respect to the solar system cannot be measured or inferred from 
pulsar-timing observations. Consequently we treat it here as a free parameter (V r > means that the pulsar 
is moving away from us), although our calculations allow us to set limits on its magnitude. 

For a specific assumed value of the radial velocity we obtain the present motion of PSR B1913+16 in an 
inertial system co-moving with the solar system. Its velocity in the Galactic (X-Y-Z) frame is then simply 

cos b cos / \ / — sin / \ / — sin b cos / 

Vq + Vt I cos b sin/ + nid\ cos/ + Hbd \ — sin sin/ I . (5) 



sin6 / V / \ cos b 



The Galactic velocity of the Sun is the Galactic velocity of the local standard of rest (see Appendix A) plus 
the peculiar velocity of the Sun (Dchnen & Binney 1997; Bienaymi]l999): 

Okm/s \ / 10 km/s 
V© = I 221 km/s + 5 km/s I . (6) 



Okm/s / \ 7 km/s 



Given the present position and velocity X p and V p , we calculate the pulsar's Galactic trajectory back 
in time using a model for the gravitational potential of the Galaxy (Appendix A). For V r < —500 km/s or 
V r > 200 km/s, we find only one intersection with the Galactic plane [Z = pc) just a few Myrs in the past 
(case A). For smaller radial velocities one finds several additional intersections, a second one ;> 35 Myrs in 
the past (case B) and even a third one ;> 135 Myrs in the past (case C). These kinematic ages, Tki n , as a 
function of V r are shown in Fig. [l]. The third intersection is well beyond the age inferred from the spin-up 
limit, r su , and the characteristic age of the pulsar, r c , and therefore we exclude this case from our analysis. 

Given our assumption that the progenitor was in the Galactic plane at the time of the second SN, we 
take the position and velocity of the intersections to correspond to the post-SN conditions. In what follows 
we consider two separate cases. In the first case (A), the binary pulsar has intersected the Galactic plane only 
once, i.e. we are seeing the direct runaway of the binary pulsar out of the Galactic plane, caused primarily 
by the kick imparted to the system at the last SN. In this case its kinematic age is ~ 3 Myr for a wide range 
of radial velocities (see Fig. |l|). In the second case (B), the binary pulsar has moved from below to above 
the Galactic plane. Note that in this case, given the oscillatory nature of the motion perpendicular to the 
Galactic plane, the time spent by the binary system close to the Galactic plane (\z\ <; 300 pc) at its present 
position, is short compared to its kinematic age (< 8%; see Fig. ||). Consequently, we consider case A as 
more likely. Nevertheless, we will investigate both cases. 



and fig have a normalized covariance of +0.3147. 



- 5 - 



If the binary system is indeed old, i.e. case B, then there is a lower limit to its age of ~ 35 Myr (c.f. 
Fig. 0). The minimum age of 60 Myr derived by Arzoumanian, Cordes, fc Wasscrman (1998) is a result of 



their assumption of constant Galactocentric radius for the motion of the B1913+16 system, which is expected 
to brake down for these long ages. 

We obtain the magnitude and direction of the center-of-mass velocity, Vcm, imparted to the system 
due to the second SN explosion by subtracting the Galactic rotation velocity from the velocity of the binary 
system at Z = pc (under the assumption that any pre-SN peculiar velocity relative to Galactic rotation 
was small). The magnitude of Vcm is given in Fig. || as a function of the radial velocity V r ., for both cases 
A and B. In what follows we use this magnitude as a constraint on the properties of the SN kick and of the 
pre-SN binary system. 

It is worth noting that none of the calculated pulsar trajectories in the Galaxy penetrates far into 
the central region, where the axisymmetric model for the gravitational potential might fail to be a good 
approximation. 



3. Orbital Dynamics at the Supernova Explosion 



In double NS systems the compactness of the two binary members excludes any mass transfer or tidal 
interaction between them after the second SN. Any evolution of the binary characteristics is due only to 
the emission of gravitational waves. For a given radial velocity, V r , we calculate the time, Tsn = Tan, 
elapsed since the second SN, i.e. the kinematic age of the binary system (see Fig. ||) and use the (first-order) 
equations of Peters |(1964) to calcul ate the post-SN semi-major axis, Aj, and eccentricity, e/. As discussed 
in more detail in Fryer fc Kalogera (1997), the various evolutionary sequences discussed in the literature for 
the formation of double NS converge to the same configuration before the second SN: a binary consisting of a 
NS and a helium star in a circular orbit (because of an earlier common-envelope or mass transfer phase that 
led to the recycling of the pulsar). We analyze the effect of the SN explosion to the orbital dynamics of the 
system and we relate the post-SN parameters (masses, orbital size and eccentricity, center-of-mass velocity, 
misalignment angle between spin and orbital angular momentum) to the pre-SN parameters (helium-star 
mass and orbital size) as well as to the magnitude and direction of the SN kick. Effects of the expanding SN 



shell on the compact companion are absolutely negligible (c.f. calculations in Fryxell fc Arnett (1981) and 



Yamaoka, Shigeyama, & Nomoto (1993), where the effective cross section of the NS companion is given by 



its radius of Bondi-Hoyle accretion). 

The (circular) pre-SN PSR B1913+16 binary system is characterized by the mass of the exploding star, 
Mi, the mass of the visible pulsar, M p , and their separation, A4. The relative orbital velocity, Vo, between 
Mi and M„ 



Vo = 



G(Mj + M p ) 

A, 



(7) 



where G denotes Newton's gravitational constant. The post-SN binary system is characterized by the mass 
of the NS born in the SN explosion, Mf = M c , the mass of the recycled pulsar, M p , the semi-major axis 
of the binary orbit, Af, its orbital eccentricity, e/, and the angle 9 between the orbital angular momenta of 
pre- and post-SN orbits. The parameters of the pre- and post-SN binary systems determine the kick velocity 



Vif imparted to the NS at the moment of its birth in the SN explosion. Using equations in Hills (1983) and 
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Kalogera (1996) we find 



V K 

X 


= ±77 




= 77 cos 




= 77 sin 



( 1 + e f)J \a{l-e f ) 



- 1 



1/2 



77 = ^ja(3(l-e 2 s ) , 



(8) 



where we used following definitions: 



Mc + M p 
Mi + M v 



(9) 



During the SN explosion, the post-SN system as a whole receives a velocity relative to the center of 
mass of the pre-SN binary. The magnitude of this center-of-mass velocity Vcm (vcm = Vcm/Vq) is given 



in |Kalogera | (1996) 



1 



< \ ,\i = «i +K2 ( 2 ) - Ks\ja{l - e 2 f ) cos9 , 



where 



M? 



M„ 2 



«1 



K 2 



K 3 = 2^KlK2 



(10) 
(11) 



{Mi + M p f ' ' (Af c + M p ) (Mj + M p ) 
If k denotes the the unit vector in direction of the orbital angular momentum of the post-SN system, then 

v cm = v cm • k = ^/kT sin 9 . (12) 

For the angle between vcm and the orbital angular momentum of the post-SN orbit we find 

a v = arccos(u|l M /ucM). (13) 

The calculations described in section 2 provide the angle between the line-of-sight, Ko, and vqm: 

i v = arccos(K • v C m)- (14) 

The angles a v , i v , and i form a spherical triangle and thus satisfy the condition 

(15) 



^ ^ COS I — COS l v COS a v < ^ 



sm i v sin a v 

Observations have determined the masses of the post-SN system, M p and M c , the semi-major axis of 
the present orbit, A, and its eccentricity e (see Table |l|). We follow back in time the evolution of the orbit 
as determined by gravitational wave emission and derive the semi- major axis, Af, and eccentricity, e/, just 
after the second SN explosion. Finally, Kramer | ( 1998) and Taylor (1999) have inferred the inclination of the 
pulsar spin with respect to the orbital angular momentum, O. Based on our understanding of the recycling 
of the pulsar, the spin axis of the pulsar (first NS) and the pre-SN orbital angular momentum axis are 
expected to have been aligned. Therefore, this observationally inferred misalignment angle, O, is in fact the 
tilt angle, 9, of the post-SN orbital plane relative to that before the explosion. 

We use all this information to constrain the (circular) progenitor system (semi-major axis A,-) and the 
birth kick velocity ~Vk (both magnitude and direction) imparted to the NS during the SN explosion of the 
immediate progenitor with mass Mj. In order to do so, we combine the values of Af, e/, and 9 with the 
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results of Section g, where we discussed the present position and (transverse) motion of the binary system 
within our Galaxy to calculate the center-of-mass velocity Vcm ; imparted to the system during the explosion 
as a function of the unknown radial velocity. Using Eq. ( |l0|) this leads to a curve in the A^-Mi parameter 
space, which is constrained by 

' < a < — — (16) 



1 + e / 1 - e / 

and Eq. 15. Based on numerical simulations of rapid accretion onto NSs ( Fryer, Benz, fc Herant~| l996), 



Fryer & Kalogera (1997) have shown that in the pre-SN progenitor of PSR B1913+16 the orbit must be 



wide enough so that the pulsar does not enter a the envelope of the helium star, otherwise it would have 
collapsed into a black hole. Therefore, there is an additional constraint imposed on (A,-,M,-) determined by 
the maximum radial extent reached by the helium stars during their lifetime. Based on stellar evolution 



calculations ( Habets 1985; Woosley, Langer, fc Weaver 1995), the maximum radius of helium stars can be 
approximated by ( Kalogera fc Webbink 1996): 



log i?Ho 



3.0965 - 2.013 log(M Ho /M ) , 
0.0557 [log(M He /M ) - 0.172]" 



2.5 



Mue < 2.5M Q 
M Hc > 2.5M Q 



(17) 



Since the maximum radius of the immediate progenitor star has to lie within the Roche lobe of the progenitor 
system, a lower limit is set on A4 for a given M; (see figure 4 in Fryer fc Kalogera~| l997), which we take into 
account. 



4. Results and Discussion 

For every allowed pair (Ai, Mi) we can now calculate the kick velocity Vjf required to form the post-SN 
system, using Eqs. ||. We are interested in its magnitude and its direction with respect to the pre-SN orbit. 
We denote the polar angle between the pre-SN orbital angular momentum and the kick velocity with ipg, 
and the azimuthal angle with ip^. (The direction ip$ = tf)^ = points away from the companion's position 
at the moment of the SN explosion.) Hence, for a given V r we get an upper and a lower limit on Ai, Mi, 
Vk, ipe, and Vv- A solution for these pre-SN parameters cannot be found for all values of V r , since values 
outside the allowed range lead to pairs (Ai,Mi) that fail to satisfy conditions (|l5|), (filf), and (fl7|). In this 
way, we are able to constrain the unknown magnitude of V r in addition to the five pre-SN binary and kick 
characteristics. 

The most uncertain parameter in our calculations is the distance to the binary system and thus its 
present location in the Galaxy. To account for this uncertainty we performed our calculations for various 
distances ranging from 5.9 to 9.7 kpc and plotted the most conservative constraints (Figs. [I] to |^). Figures 
^ and U correspond to case A, where the pulsar is young (~ 3 Myr), with 9 — 22° (prograde) and 9 = 158° 
(retrograde), respectively. Figures || and ^ give results for case B, again for the prograde and retrograde 
cases. Both the mass, Mi, of the exploding helium star and the pre-SN orbital separation, Ai, are well 
constrained in the retrograde cases A and B and the prograde case B. The kick magnitude, Vk, is better 
constrained in case B, while both retrograde cases require higher kick magnitudes, as expected. For all 
four cases the direction of the kick is narrowly constrained to be very close to the pre-SN orbital plane and 
opposite to the orbital motion of the exploding helium star. In the retrograde cases the angles and ips 
are restricted to extremely narrow ranges. 



According to Kramer 1998 the observational data for the pulse profile evolution are equally well described 
by a prograde and retrograde system. He however argued that a retrograde system would require typically 
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larger and therefore less probable kick velocities than for the prograde solution, and therefore considered 
this case as less likely. This tendency is indeed visible in the results of our calculations (Figs. || to[7|). More 
importantly, a comparison between the solid angle covered by the direction (ipei^^) m the prograde and 
retrograde cases (for each allowed value of V r ) indicates that the retrograde cases require a much higher- 
degree of fine-tuning in the kick direction. For this reason we consider the retrograde cases to be less 
probable. 

It is important to point out that in all four possible cases the direction of the kick imparted to the newly- 
born NS is strongly restricted to be almost in the orbital plane of the progenitor system (ipg w 75° — 85°). 
Since the spin of the helium star is expected to be aligned with the pre-SN orbital angular momentum due to 
the mass transfer in the binary, our results imply that the kick was given in a direction almost perpendicular 
to the spin axis of the NS progenitor. Such a constraint has interesting and restrictive implications for the 
physical mechanism responsible for the kick in terms of the relevant direction and time scale of the process. 
Whatever the details of the mechanism, it must be such that it can allow the kick to be directed almost 
perpendicular to the rotation axis of the collapsing star. Further the time scale relevant to the kick mechanism 
must be short compared to the birth spin period of the NS. Otherwise kick components perpendicular to 
the NS rotation would be cancelled resulting in a kick parallel to the NS spin axis, which contradicts the 



constraints derived here. Our findings are in agreement with Birkcl & Toldra 1997, Cordes & Chernoff 1998 



Deshpande, Ramachandran, & Radhakrishnan 1999 who argue that there is no correlation between the pulsar 
spin and velocity vectors. We also note that if the kick mechanism is related to the escaping neutrinos during 
the collapse, then the spin period of the proto-NS must be longer than the neutrino-diffusion timescale (a 



few seconds, Burrows fe Hayes 1996) 



The radii of helium stars in late stages of their evolution are not well defined quantitatively, although 
the qualitative result of rapid expansion of relatively low-mass helium stars is common to all the published 



calculations (Habets 1985; Woosley, Langer, fe Weaver 1995). We have explored the sensitivity of our results 



by artificially decreasing i?He,max by a factor of 2. Quantities describing the kick magnitude (vk ) and the 
kick direction (tpg, ip^) turned out to be rather insensitive to changes of i?Hc,max- However, the lower bounds 
on the orbital separation, Ai, and the mass of the Helium star, M^, just before the second SN explosion 
proved to change significantly if one reduces i?He,max by a factor of 2. The lower bound on Ai decreases by 
-50% for the 9 = 22° cases (Figs. |§) and -20% for the 9 = 158° cases (Figs. |,@). The lower bound on 
Mi decreases by —20% for the 9 — 22° cases, while for the 9 — 158° cases the changes in the lower bound 
on Mi are less than 15%. 

In section || we calculated the kinematic history of the PSR B1913+16 system based on the model for 



the Galactic gravitational field of Kuijken fe Gilmore (1989). In particular, calculations done for the case B 



( T kin ^> 35 Myr) could be sensitive to deviations of the model from the true Galactic potential. Therefore 



we repeated all our calculations using the Galactic potential model published by |Paczyhski | (1990). The 
resulting changes in kick and orbital parameters proved to be below the 10% level. We further tested our 
assumption of the PSR B1913+16 system being born at Z = pc, using in addition Z = ±100 pc as the 
Galactic height at birth. Our analysis turned out to be adequately robust against these changes. 

Finally, we can use the upper limit on \V r \ to constrain the uncertainty in the NS-mass measurements. 
As pointed out by Damour fc Dcrucllc | (1986), timing observations of binary pulsars do not allow the 



determination of the intrinsic masses of the binary system. Instead the (intrinsic) masses divided by the 
Doppler factor, which to first order is given by D ~ 1 + V r /c, are measured. Thus our limit of \V r \ < 
1200 km s" 1 implies \M ohs /M - 1| < 0.004. This limit on the mass uncertainty could be improved if an 
upper bound on the helium-star mass leading to NS formation is obtained theoretically (from SN collapse 
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calculations) in the future (see Fig. ||). 

Although some of the limits imposed on the parameters of the immediate progenitor system and the 
birth kick are still rather wide, improved values for the distance to PSR B1913+16 would certainly tighten 
the constraints. In the future an improved distance estimate may be derived from refined measurements of 
the relativistic orbital-period decrease (c.f. Bell & Bailcs 1996). Also, a better estimate for the age of the 
binary system or a reliable mass limit on the helium progenitor star from theoretical considerations combined 
with the calculations here, would lead to clearly stricter limits. 
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A. The Galactic potential 



It is convenient to use standard cylindrical coordinates (R, Z, <p) with the origin at the Galactic center, 
R is the Galactocentric radius, 4> is the azimuthal angle and Z is the height above the Galactic plane. The 
potential considered here is axi-symmetric about the Z-axis and symmetric with respect to the Galactic 
plane, i.e. $ = <£>(i?, \Z\). Following |Carlberg fc Innanen~| (1987) and [Kuijken fc Gilmore] (1989) the Galactic 
potential is a sum of three components: 



where 



&a(R, Z) 



'I'Gal = $1 + $2 + $3 
GM A 



A = 1,2,3 



(Al) 
(A2) 



The numerical values for the parameters Ma, a, b, (3f, and hf are given in Table |[ The equations of motion 
for a star in the Galactic potential $ = <I>Gai {R, Z) arc 



■■ _ X <9$ ■■ _ Y d<f> •• _ 

A ^ ^ Y -RdR' Z ~dZ 



(A3) 



R dR ' 

In the equations of motion we changed to Cartesian coordinates in order to avoid the singularity of cylindrical 



coordinates at R — 0. The equations of motion (A3) are then solved numerically using a fourth-order Runge- 



Kutta routine (Press et al. 1986). 

The speed for a circular orbit in the Galactic plane at a distance i?» from the Galactic center is 



V rot (R* 



5$ 



(A4) 



At the solar distance, i?* = 8.5 kpc, the model here gives a circular speed of 221 km/s for the local standard 
of rest. 
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Fig. 1. — Kinematic age, Tki n , of the PSR B1913+16 system as a function of the (unknown) radial velocity, 
V r , for the three cases A, B, and C. The characteristic age of the pulsar, t c , is indicated by a dashed line. 
The dotted line indicates the age, r su , inferred from the spin-up limit during the accretion phase (taken from 



Arzoumanian, Cordes, fc Wasserman 1998) 



Fig. 2. — Time At required for the PSR B1913+16 system in case B to move from z — —300 pc to z — +300 
pc devided by its kinematic age as a function of radial velocity V r . The time spent close to the Galactic 
plane is less than 8% of the kinematic age of the PSR B1913+16 system. 

Fig. 3. — Velocity of the PSR B1913+16 system with respect to its Galactic standard of rest at the moment 
of birth [Z — pc) as a function of the (unknown) radial velocity, V r , for case A (dashed-dotted line) and 
case B (dashed line). 

Fig. 4. — Limits on the immediate progenitor (helium-star) mass, Mj, the pre-SN orbital separation, A4, 
the magnitude, Vk, of the kick velocity imparted to the pulsar companion during the SN explosion, and its 
direction given by the two angles, ip<p and ipg, as a function of the (unknown) radial velocity, V r , for case A 
and 9 = 22°. The largest regions (most conservative results) were obtained for the maximum distance d of 
9.7 kpc. 

Fig. 5. — Same as in Fig. || for case A and = 158°. The largest regions (most conservative results) were 
obtained for the maximum distance d of 9.7 kpc. 

Fig. 6. — Same as in Fig. ^ for case B and 6 — 22°. The largest regions (most conservative results) were 
obtained for the maximum distance d of 9.7 kpc. 

Fig. 7. — Same as in Fig. [| for case B and 9 — 158°. The largest regions (most conservative results) were 
obtained for the minimum distance d of 6.9 kpc. 
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Table 1: Parameters for the PSR B1913+16 system 



R. A. (J2000) 


a 




19 h 15 m 28?00 


Decl. (J2000) 


S 




16°06'27'/4 


Proper motion in R. A. (mas yr _1 ) 


Ha -- 


- a cos 5 


-3.27 ±0.35 


Proper motion in Decl. (mas yr _1 ) 


M<5 = 


-- 5 


+1.04 ±0.42 


Distance (kpc) 


d 




8.3 ± 1.4 


Spin period (ms) 


P 




59.03 


Period derivative 


P 




8.63 x 10~ 18 


Characteristic age (Myr) 


T c = 


: P/2P 


109 


Mass of the pulsar (M©) 


M p 




1.44 


Mass of the companion (M©) 


M c 




1.39 


Present semi-major axis (i?©) 


A 




2.80 


Present eccentricity 


e 




0.617 


Orbital inclination (deg) 


i 




47.2 
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Table 2: Parameters defining the Galactic potential 







M A 
(1O 9 M ) 


a 

(kpc) 


b 

(kpc) 


at 


hi 
(kpc) 




hi 
(kpc) 




hf 

(kpc) 


A = 


1 


145 


2.4 


5.5 


0.4 


0.325 


0.5 


0.090 


0.1 


0.125 


A = 


2 


9.3 





0.25 




















A = 


3 


10 





1.5 






















V r (km/a) 



